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Abstract. We calculate the radiative decay widths of the ao(980) and /o(980) scalar mesons into p7 and 
luj considering the dynamically generated nature of these scalar resonances within the realm of the Chiral 
Unitary Approach. The main ingredient in the evaluation of the radiative width of the scalar mesons are 
the loops coming from the decay into their constituent pseudoscalar-pseudoscalar components and the 
subsequent radiation of the photon. The dominant diagrams with only pseudoscalar mesons in the loops 
are found to be convergent while the divergence of those with a vector meson in the loop are written in 
terms of the two meson loop function easily regularizable. We provide results for all the possible charge 
channels and obtain results, with uncertainties, which differ significantly from quark loops models and 
some version of vector meson dominance. 



PACS. 14.40. Cs - 13.25.Jx - 13.40. Hq 

1 Introduction 

The radiative decay of resonances has been advocated 
as a privileged tool to study the nature of resonances 
[HE]- A good example of the interest in these reactions 
is offered by the large amount of works devoted to the 
study of the radiative decay of the Z?* (2317) charmed 
scalar meson [3j. The uncharmed scalar mesons are the 
object of much debate regarding its nature as a possible 
qq state, a meson-meson molecule, a dynamically gener- 
ated object from the interaction of pseudoscalar meson 
pairs in coupled channels or a hybrid of those structures 
[4 l l5 ] l6 l [7]. The idea that the light scalar mesons cr(600), 
/o(980), a (980), k(800) are bound states or resonances 
formed from the interaction of pseudoscalar mesons has 
been gaining support, particularly because the underlying 
interaction is well known from the chiral Lagrangians [5], 
and then using any reasonable nonperturbative tools, like 
the Bethe Salpeter equation [T Wl llO p ilj. the N/D method 
[12] or the Inverse Amplitude Method 13,14 the scalar 
states appear automatically. The support for this picture 
is also strengthened by the success reproducing different 
reactions where these resonances are produced without 
the need to introduce any free parameters. Examples of 
this are the description of the <j> — ► 7r°7r°7 reaction with 
its strong /o(980) peak, or the <f> — > tt°7]j reaction with 
its prominent do (980) peak [1511101116). the 77 — > -k-k and 
77 — > KK reactions [TTlfTB] and, adjusting some normal- 
ization parameters, the description of the J '/& — > NNtttt 



[19] and the J/W — > ujtttt, J/ty — ► <fm-K, with clear sig- 
nals for the cr(600) in the first reaction and of /o(980) in 
the second one [20II2T] . The apparent narrowness of the 
er(600) in the J/W — > uitit reaction is nicely interpreted 
in since the strength of a production is proportional 
to the ratio t/V, with t the scalar isoscalar ww amplitude 
and V the potential in the same channel (the Adler term) . 
The narrowness is due to the strong energy dependence 
of the Adler term, which grows faster that the scatter- 
ing t-matrix. Further considerations concerning the scalar 
meson sector are done in [22) . 

The success of the chiral unitary approach in describ- 
ing the scalar sector and in particular the /o(980) and 
ao(980) resonances, can be further tested with new ob- 
servables like those we suggest in the present work. Given 
the fact that the <f> — > 7r°7r°7 and 4> — ► ir i]"f reactions 
are well understood within this picture, one can look at 
the crossed channel reactions of the former ones in places 
where the energy conservation allows it. This is the case 
of the /o(980) and cio(980) resonances decaying into Vj 
where V is now a vector meson which replaces the <fi 
to make the reaction energetically possible. Such vector 
mesons can only be the p and the u>. The sensitivity of 
this observable to details of different models is huge and 
has been discussed in several works, assuming vector me- 
son dominance and direct coupling of the scalar to VV 
[23j . or quark loop contribution 1 or meson- meson loop 
contribution, KK Ij, and KK and 7T7t [21] • The work 
presented here will follow the ideas of the chiral unitary 
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approach used to describe the <fi — > 7r°7r°7 decay, evalu- 
ating loops of KK and tttt, as in the former works, but 
also those containing one vector meson, which proved to 
be relevant in the study of the radiative decay of vector 
mesons to two pseudoscalar mesons and one photon J35J 
[T6] . An experimental proposal to measure the /o, floj to 
V-y decays has been approved at COSY/WASA [26]. 

2 Formalism 

By using the techniques of the unitary extensions of chi- 
ral perturbation theory (chiral unitary approach or Unita- 
rized ChPT) the low lying scalar resonances (a, k, /o(980) 
and ao(980)), among many other mesonic and baryonic 
resonances, are generated dynamically. With the only in- 
put of the low-lying chiral Lagrangians, the implementa- 
tion of unitarity in coupled channels and the exploitation 
of the analytic properties of the scattering amplitudes, the 
approach generates poles in unphysical Riemann sheets of 
the unitarized meson-meson scattering amplitudes which 
can be associated to those resonances. Hence they qual- 
ify as quasibound states of their constituent meson-meson 
components. Furthermore, by evaluating the residue of the 
meson-meson scattering amplitudes at the pole positions, 
this picture provides the value, including the phase, of the 
couplings of these resonances to their constituent building 
blocks. Indeed, close to the pole position the Laurent ex- 
pansion of the scattering amplitude in a particular isospin 
and partial wave can be approximated by its dominant 
term 



where i and j refer to a given meson-meson estate and 
s p ~ (M — iT/2) 2 is the pole position, with M and r the 
mass and width of the associated resonance. Hence, gi in 
Eq. ([1]) can be identified as the coupling of the dynam- 
ically generated resonance (with the quantum numbers 
of the amplitude t) to the i — th channel. In particular, 
using the Bethe-Salpeter unitarization procedure, we ob- 
tain the couplings of the ao and fo to their constituent 
pseudoscalar-pseudoscalar (PP) components shown in ta- 
ble [T] The difference with the results of ref. [27] will be 
considered as an uncertainty in the error analysis in the 
present work. 
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\ 9i \ from ref. 27 


f -mi 


469 - 


il247 


1332 


1100 


f KK 


-3369 - 


il606 


3732 


3680 


a KK 


-4732 


-i943 


4825 


5500 


a ( )7T?7 


3166 


- i712 


3245 


3900 



Table 1. Couplings of the /o(980) and a (980) to their dif- 
ferent pseudoscalar-pseudoscalar constituent channels. All the 
units are in MeV. 



Once we have these couplings, the philosophy we fol- 
low in order to obtain the radiative decay widths is al- 



most straightforward: one has to consider the transition 
from the scalar mesons to the possible PP at one loop and 
then attach the photon to the possible allowed places, con- 
sidering that a vector meson in the final state needs to be 
produced. For reasons that will be clear later on, the di- 
agrams that we need to evaluate are only those shown in 
ng.ffl 

The channels we will consider are fo — > p 7, fo — > u>j, 
clq — ► p°7, <2q — > W7 and — > p + 7. In fig. [l] S represents 
the decaying scalar mesons; Pi, P2, pseudoscalar mesons 
and Vi, V2, vector mesons. In the Appendix, tables[3][71 we 
show the different allowed P1P2V1V2 particles of the dia- 
grams in fig.[T]together with the corresponding coefficients 
for each channel to be explained latter on. In fig.[l] P, q, k 
and Q represent the momentum of the different lines that 
will be used in the evaluation of the loop function. 

At this point we take advantage of the implications of 
gauge invariance to simplify the calculations. We follow 
a similar procedure as done in refs. |28U15j in the evalu- 
ation of the radiative <fi decay and in refs. [291130] for the 
radiative axial- vector meson decays. 

The general expression of the amplitude for the radia- 
tive decay of a scalar meson into a vector meson and a 
photon (S — ► V7) can be written as 

T = e v ^T^ (2) 

with 

T^ v = a g» v + bQ^Q u + cQ»k" + dk^Q u + ek^k u (3) 

where Q is the final vector meson momentum and k the 
photon momentum, which are the only independent avail- 
able momenta. In Eq. ([2]), ty and e are the vector meson 
and photon polarization vectors respectively. On the other 
hand, due to the Lorenz condition, ey^Q^ = 0, t v k v = 0, 
all the terms in Eq. ([3]) vanish except for the a and d terms. 
Furthermore, gauge invariance implies that T^ v h v = 0, 
from where one gets 

a = -dQ-k. (4) 

Therefore the amplitude gets the general form 

T = -d(Q.k gT - k»Q»)e Vll e v . (5) 

Hence, we only need to evaluate those diagrams contribut- 
ing to the d-teim, that is, those having a final structure 
k^Q v ' . The advantage to evaluate only the d coefficient 
is that only the loop diagrams of fig. [1] contribute since 
other diagrams, like those involving photon couplings to 
the vertices which are necessary to fulfill gauge invari- 
ance, do not give contribution to the d coefficient [28,15, 
[29] . Another advantage is that, from dimensional reasons 
(performing explicitly the Feynman integrals), we will see 
that the d coefficients are finite for the type-a diagrams 
of fig. [TJ while the logarithmic divergence from the type-b 
diagram can be easily identified and regularized by ex- 
pressing it in terms of the two meson loop function of the 
PP scattering problem. 

With the previous discussion in mind, let us explain 
the explicit steps in the calculation of the radiative decay 
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(Q-q) 

a) b) 

Fig. 1. Feynman diagrams contributing to the scalar radiative decays 



width. We will consider first the evaluation of the type-a 
diagram of fig. [U The Lagrangians needed in the evalu- 
ation of the type-a diagram in fig. [T] are, for the vector- 
pseudoscalar-pseudoscalar (VPP) vertex: 



: VPP = -iJL{v>[B l j>,p\; 



(6) 



where g — —MyGy/f 2 , with My the vector meson mass 
and Gy a coupling constant defined in [31) and for which 
we use the numerical value Gy = 55±5 MeV from ref. [16] . 
In Eq. ([6|), V and P are the SU(3) matrices containing 
the octets of vector and pseudoscalar mesons respectively, 
see e.g. ref. [31] . and () means that one has to perform 
the trace of the SU(3) resulting matrix. In Eq. ©, / is 
the pion decay constant (/ = 93 MeV), however one can 
assign an uncertainty to the / constant since it could range 
from fa to fr). For the calculations we will actually use, as 
central value, / = 1.08 x 93 MeV and we will also consider 
uncertainties in our calculations by taking a range of / 
between /„■ and 1.15/,,. 

The PP7 vertex can be readily obtained from the low- 
est order ChPT Lagrangian [8] 



from where the PP7 amplitude follows: 



it 



PP 7 



(7) 



(8) 



with pi(j)2) the momentum of the ingoing(outgoing) pseu- 
doscalar meson, Qi is the sign of the charge of the pseu- 
doscalar meson and e is taken positive. 

Therefore, the amplitude for the type-a diagrams takes 
the form: 



it„ 



-iAg S p 1 p 2 



d q i 
(2tt) 4 {q + kf -mj + ie 

i 



q 2 — m\ + ie (Q — q) 2 — m\ + ie 
x {-i)eQ x e v {q + q + kf 

x (-i) C ^fe v ,{ q -Q + qr, 



(9) 



where A are coefficients given in table [5] needed to relate 
the SP1P2 coupling in charge basis with those in isospin 



basis and with the unitary normalization of refs. [7j, gsp t p 2 ■ 
In Eq. (jQj) Qi is the sign of the charge of the Pi pseu- 
doscalar meson, 7711(7712) is the mass of the PiiP^) pseu- 
doscalar meson, My the mass of the final vector meson 
and C are coefficients coming from the Lagrangian of 
Eq. ([6]) after performing the trace of the matrix and which 
depend on the particular Pi, P2 and V particles (specifi- 
cally, C is the coefficient of (V^P, P]) of Eq. © defined 
as C(P 1 dP2-P 2 dP 1 )). The different A, C, coefficients are 
given in the Appendix, table [5] 

By using the Lorenz condition, ey^Q^ = and t v k v — 
0, the amplitude reads 



^ „M v Gy „ 
to = iAg SPl p 2 eQ 1 C ^_ iey^ 



d 4 q 



X / t^WV 



1 



1 



(27r) 4 "* ^ (q + k) 2 -ml+ ie q 2 - m\ + ie 
1 



{Q - q) 2 - m\ + ie 



(10) 



This loop integral can be easily done by using the Feyn- 
man parametrization. We use the identity 



1 _ r 1 r x 1 

abc J Q J [a + (b - a)x + (c - b)y] 3 ' 



(11) 



We set 

a = (Q - q) 2 - ml, 

u 2 2 

= q — m 1 , 

c = (q + k) 2 — m\ 

and perform the change of variable 

q = q' + (1 - x)Q - yk. 



(12) 



(13) 



Now we have to recall the above discussion in order to keep 
in mind that we only need the d coefficient (see Eqs. (O 
and ([5])). That means that we only need the terms pro- 
ducing a final structure of the type k^Q" , which reduce 
considerably the number of terms and, more important, 
the resulting expression contributing to the ci-term is just 
finite. Thus the apparent superficial logarithmically di- 
vergent type-a loop (see Eq. ([9]) becomes pure convergent 
and, hence, with no need for rcgularization, it is univo- 
quely defined (no regularization parameters, scale, etc). 
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For the type-a loop, the final, simple and finite, expres- 
sion for the d-coefficient of Eqs. ([3]) and ([5]), from where 
the amplitude can be obtained by using Eq. ([5]), is 



d a = AC Qi g S p 1 p 2 V2e 



MyGy 1 
P 8tt 2 



dx / dy 



(1 ~ x)y 
s + ie 
( 14 ) 

where s = Q 2 x(l — x) + 2Q- fc(l — x)y+ (m| — m\)x — m 2 . 

In the final step of the derivation of Eq. (|T4]) . since the 
remaining d 4 q' integration is finite, we have used that 32J 



- ie) 3 



is 



(15) 



where G' = 3.g' 2 /(47r 2 /) with g' = -GyM p /{V2f 2 ). Since 
in the type-b loops we have two vertices of the type VVP, 
the amplitude is proportional to G' 2 or g' A . Hence, the 
contributions to the decay width of the type-b loops go 
like g' s . This means that small differences in the value of 
g' are strongly magnified in the evaluation of the radia- 
tive decays widths from the type-b mechanisms. Therefore 
a good numerical value for the VVP couplings is called 
for. Hence, in order to fine tune the numerical value of 
the VVP coupling constant we proceed as follows. From 



So far, what we have done corresponds to the me- 
son loop calculation of [T|l24] with a significant difference, 
which is the coupling of the scalar resonance to the two 
pseudoscalar mesons. In [I] it is taken from [33] obtained 
using the Weinberg compositness condition 34J, and as 
quoted in [Tll2"4"] the estimation of the coupling should be 
viewed as qualitative [sic] pQ. Furthermore, only KK loops 
are contained in pQ. In 7T7t loops are also considered 
and the couplings of the scalar to the pseudoscalar mesons 
are taken from the effective resonance Lagrangians of [31] . 
One should note in this respect that, while the vectors can 
be considered as genuine new fields, additional to the pseu- 
doscalar ones, the scalar fields cannot be considered at the 
same level since they come from the same pseudoscalar- 
pseudoscalar Lagrangians after a proper unitarization 7, 
[T2 1 I91IT0 1 ITT] . This means that the use of a proper unitary 
theory involving the pscudoscalar-pseudoscalar meson chi- 
ral Lagrangians and the scalar resonance Lagrangians of 
[3T] leads to double counting. On the other hand, assuming 
that these scalar resonance Lagrangians are only a means 
of providing at tree level the couplings of the scalar to 
pseudoscalar mesons, there are also problems since the ao 
and fo properties cannot be simultaneously fitted with 
the structure of the resonance Lagrangians of [31] . In our 
case, as explained above, the couplings are taken from 
the residues at the poles of the unitarized meson-meson 
interaction amplitudes, which reproduce very accurately 
the experimental data in a range of energies from thresh- 
old to 1.2 GeV and are consistent with a large variety of 
physical processes. 

The other novelty of our approach is the evaluation of 
the type-b loops. The idea to include these loops stems 
from the relevance of intermediate vector mesons in the 
radiative decay of p, u>, into meson-meson-photon |25| . 

For the evaluation of the type-b loops, we first need 
the V"f vertex given by 



ty-y — —FyeXyMy ty ■ e 



7> 



(16) 



where Ay is 1, 1/3, — \/2/3 for V — p, v, <fi respectively, 
Fy is the coupling constant in the normalization of [31] 
for which we use the value Fy — 156 ± 5 MeV [15] . 
For the VVP vertex we use the Lagrangian [351136] : 

C VVP = ^=e^ af) (d^d a V p P), (17) 



V V 

Fig. 2. V — * Py decay diagram 



Eqs. IT6]) and ([17]) . the decay width of the decay of a vector 
meson into a pseudoscalar meson and a photon, V — > P7, 
(see fig. E]), takes the form: 



e 

247 



Bj G'FyX\ 

My. 



(18) 



where F>i are numerical coefficients depending on the dif- 
ferent VVP channels that will be explicitly defined later 
on, and |fe| is the center of mass momentum of the photon 
or the pseudoscalar meson. In table[2]we show the theoret- 
ical branching ratios together with the experimental ones 
[37] . The agreement with the experimental data is fair but 



i 


BRf 


BR e X p 


p ± TY ±r y 


3.6 x 1(T 4 


(4.5 ±0.5) x 10~ 4 




3.6 x 1(T 4 


(6.0 ±0.8) x 10" 4 


p u vi 


3.1 x 10" 4 


(2.7 ±0.4) x 10~ 4 


u 

con 7 


6.06% 


8.91 ± 0.24% 


w?77 


6.5 x 1(T 4 


(4.8 ±0.4) x 10~ 4 




9.3 x 10" 4 


(9.9 ±0.9) x 10~ 4 




19.1 x 10" 4 


(23.1 ± 2.0) x 10" 4 



Table 2. Theoretical and experimental branching ratios for 
different vector meson decay processes. 



the results can be improved by incorporating SU (3) break- 
ing mechanisms. For that purpose, we will normalize here 
the G' coupling by multiplying it by N l = y/BR^ xp /BRf, 
dependent on the particular channel, such that the theo- 
retical branching ratios agree with experiment. For the 
evaluation of the second column of tableland in the rest 
of the paper we fix G' = 0.098 MeV -1 . The experimental 
errors in table [5] will be considered later on when evalu- 
ating the uncertainties of our final results for the scalar 
radiative decays. 
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The amplitude for the type-b diagram of fig. [I] reads 



d 4 q . NpG' 



it b = -iAg SPl p 2 J j^-^iB—j=-F v M Vl e\ Vl j^ 



Vi 



V2 



(q + k) 2 — m\ + ie (Q — q) 2 - m\ + i 



x — — I -gppi 



q 2 - M 2 + it 



Ml 



(19) 



where Mi (M2) is the mass of the V\ (V2) vector meson. In 
Eq. lfl9|) . A has the same meaning as in the type-a loop 
case and B is the coefficient of the P1V1V2 vertex obtained 
after performing the SU(3) trace in (VVP) (see Eq. (]17p) 
defined as B P\V{V2- Analogously, C is the coefficient com- 
ing from the P2V2V vertex defined as C P2VV2 from the 
resulting expression after taking the trace in (VVP). The 
Nb and Nc couplings are the normalization factors for 
the B and C VVP vertices just that the V — ► Pj decays 
agree with experiment, as explained above. In the eval- 
uation of the type-b loops we use for each VVP vertex 
a value of Ni obtained from the average of the different 
charge channels of table [2] containing the same vectors. 

The loop integral in Eq. (|19j) is apparently quadrati- 
cally divergent, the highest degree of divergence coming 
from the last qpqpi term from the propagator of the vec- 
tor meson in the loop. However, this last term gives zero 
contribution since q a qa' qpqp' is a totally symmetric tensor 
and vanishes when contracting it with the antisymmetric 
Levi-Civita tensor. 

Thus, the remaining logarithmically divergent ampli- 
tude reads: 



N B G' 1 N C G' 

-iAg S p 1 p 2 B — —F v e\ Vl ^^C- 



V2 



Vi 



V2 



d A q 

WT 4 



ftvap Ai'v'a! u o 



1 



(q + k) 2 -m\ + ie (Q - q) 2 
1 



r 



M 2 



- ic 



(20) 



Even though the previous expression is logarithmically 
divergent, the divergence can be isolated and expressed 
in terms of the two-meson loop function that appears in 
the unitarization procedure of the meson-meson scattering 
amplitudes. The reasoning is as follows: by using that 



^(T p = 



n"M' „w' „ua' 



g OCU g^ 



(21) 



the e ^p e ^'u' a ' k^Q^q^ e v e Vv , term of Eq. §M) is 
reduced to 



X {k ■ Q q 2 g^ + q ■ Qk^q v + q ■ kq^Q" 

- q-kq-QgT - q 2 k"Q u -k-Q q^q") . (22) 

In order to separate the divergent part in a convenient 
way, we write the bracket {...} of Eq. (|22"]) . by adding and 
subtracting q 2 (k ■ Qg^ v — k' J -Q u )/2, in the following way: 

{...} = [A-QgV + + q-kq»Q u 

- q-kq-Qg^ - q 2 k»Q v -k-Q q»q v 



- \k-Qq 2 g^ + \q 2 ^Q v ] 
+ \q 2 (k-Qg^ -WQ V ). 



(23) 



Now the square bracket [...] term in Eq. (|2"3]) leads to a 
finite contribution after performing the Feynman parametriza- 
tion since the divergent contributions coming from the q' 
terms of Eq. (fni|) cancel algebraically. 

The remaining part coming from the last q 2 (k ■ Qg^ v — 
k fl Q u )/2 term of Eq. (|2"3"j) leads to a logarithmically di- 
vergent contribution which can be written in terms of 
the known two-meson loop function, G(s, m%, 7712), already 
used in the meson-meson unitarized scattering amplitude. 
Explicit expressions for G(s, m\, m-i) can be found, e.g., in 
refs. [381139], Indeed, we can do the following transforma- 
tion when considering the propagators present in Eq. I|20p : 



\q 2 (k-Qg^ -k»Q»)^ 2 



1 



q 2 -M 2 
1 



(q + k) 2 — m{ (Q — q) 2 — m\ 



q 2 -M 2 



MJ 



(k ■ QsT WQTj { ?2 _ ^ , q2 _ ^ 
1 



(q + k) 2 — m\ (Q — q) 2 — m\ ' 



(24) 



The term 



Mi 



-Mr, 



in the big parenthesis leads to a conver- 



gent part while the 



-Mi 



term is proportional to the two 



<? 2 -M 2 a 

meson loop function, G((Q + k) 2 ), which can be properly 
regularized either with a cutoff [7J or with dimensional reg- 
ularization [12,13]. The procedure followed here to write 
the divergent part in terms of the meson-meson loop func- 
tion is similar to what was done in ref. [40) in a different 
context. 

Therefore, gathering the convergent parts and keeping 
only the contributions to the d— term of Eq. we get 
the following convergent d-coefficient: 

^on ABCn N B N C G' 2 F V 1 
"6 = -ABC 5sPiP 2 7777 e *Vx 



2M X 



32ir 2 

2 )>/f2\ 



, 1 dx dy -^(cy(i- x y-Mj) (25 ) 

Jo Jo s + l£ 
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where now s' = Q 2 x(l — x) + 2Q-k(l — x)y J r {m\ — M$)x+ 
(Mj - ra\)y - m\. _ 

while the contribution to the d— term from the divergent 
part is 



= -ABC g SPl P % 



N B N C G' 2 F V 
4Mi 



e\ Vl G(P 2 ,mi,m 2 ) 



The total amplitude for the radiative decay process is 
then obtained from Eq. © where d = d a + d£ on + d^ lv 
from Eqs. QU), ([25]) and ([26). 

The radiative decay width of the scalar resonances into 
a vector meson and a photon is then readily obtained, 
taking the narrow resonance limit in a first step, by 



r(M s ,M v ) = 



(27) 

where Ms stands for the mass of the scalar meson. 

In order to take into account the finite width of the 
scalar resonance and the final vector meson we fold the 
previous expression with their corresponding mass distri- 
butions: 



1 



(M s +2r 3 ) 2 



ds<- 



\r 2 / " ,Ss 
- /V7r J(M a -2r g y J{M v -ir v Y 



(M v +2r v ) 2 



ds\ 



x Im 



x Im 



ss 



Ml 



iM s r s 



i 



), 



s v -M 2 + iM v r v 

' (28) 

where is the step function, s^y*. is the threshold for the 
dominant S(V) decay channel and Af is the normalization 
factor of the spectral distributions: 



TV = 



1 



(M s +2r s 



(M s -2r s ) 2 



ds<- 



x Im 



MJ + iM s r s 



(M v +2r v ) 2 
(M v -2r v ) 
Im 



ds\ 



ss 



s v 



M 2 



iM v r v 
(29) 



Advancing some results, the effect of the folding with 
the vector-meson spectral function has little influence in 
the radiative decay widths. On the contrary, the convo- 
lution with the scalar meson mass distribution is crucial 
for the decays of the /o, as we will explain in the Results 
section. 



3 Results 



tiorQ. The numbers besides the labels tttt, KK and -kt] 
are the decay widths that we would obtain had we used 
only the loops-a or -b where the scalar meson vertex is at- 
tached to 7T7T, KK and irr) respectively. The other number 
in columns 2 and 3 is the global loop-a or loop-b contri- 
bution. The theoretical errors quoted in our final results 
have been obtained by doing a Monte-Carlo sampling of 
the parameters of the model within their uncertainties as 
have appeared along the text. We have checked that the 
largest source of error in the final results for the fo decays 
is the uncertainty considered in the pseudoscalar decay 
constant /, while the largest one in the ao decays comes 
from the uncertainty in the g ao pp couplings. 

For the fo — > p°j the irir type-a loops contribution, 
in spite of being only a 10% of the KK one, influences 
the global loop-a contribution by about 30% due to the 
interference with the dominant KK type-a loop. 

In table[3]we can see that the contribution of the type- 
b loops to the decay widths, evaluated for the first time 
in the present work, is quite important for most of the 
decays, particularly for the fo — > p 0r y, ao — > fry and ao — ► 
uj"f decays. The KK loops in the type-b diagrams are small 
by themselves with respect to the irir or irr) loops. However, 
they are important to give the global type-b contribution 
after the coherent interference with the tttt or irr) loops. 
For instance, had we only considered the type-a loops, the 
results would be the same for the ao — * fry and ao — > W7 
radiative decay widths. (The differences in table [3] are due 
to the different masses of p and w). The type-b loops for 
these ao decays are dominated by the irr] loops. However, 
the irr] mechanisms of the type-b loops for the ao — > W7 
decay is one order of magnitude larger than that of the 
a o - * fn- This is essentially due to the fact that the vector 
meson attached to the photon (Vi in fig. [TJs) is a p in 
the former case and an u) in the later one. This implies a 
factor three difference in the Ay x coefficient of the direct 
coupling of the vector meson to the photon (see Eq. (fl6|) ). 
This makes, after the interference with the other pieces, 
that the width of the ao — » loj decay is much larger than 
the ao — > P7 decay width. 

Since all these interference effects are so important, 
it is crucial to know the relative sign among the different 
couplings of the model. At this point the chiral unitary ap- 
proach makes a decisive contribution since, as explained at 
the beginning of section^ it provides the gspp couplings 
including their phase. Hence the right direction of the in- 
terferences among the different mechanisms is a genuine 
non-trivial prediction of UChPT. 

On the other hand it is worth mentioning the impor- 
tance of the convolution with the scalar mass distribution, 
Eq. (|28|) . Had we not considered the folding of the decay 
width with the scalar meson spectral function we would 
have obtained radiative decay widths for the fo about 3 
times smaller than in the case when the convolution is 
performed. Indeed, as already mentioned in refs. [1 . 24J the 
loop functions for the type-a diagrams are strongly depen- 



In table [3] we show the contributions of the type-a and 
-b loops to the radiative decay widths under considera- 



1 By ao — > P7 we denote both a[] — > p°7 and 
the results are the same for both decays. 



p ± 7 since 
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loop a 


loop b 


total 


fo - P°7 


tttt: 0.44 
ifif: 4.10 


tttt: 1.82 
JOf: 0.56 


4.2 ± 1.1 


fo -> ^7 


ifi? 4.31 


™ L34 i 49 
Kff: 0.013 


4.3 ± 1.3 


a -> P7 


A'i? 7.43 


m 0.68 
TfA": 0.046 


11 ±4 


ao — » t^7 


KK 7.85 


7^:4.69 
KK": 0.79 


31 ± 13 



Table 3. Contribution of the different mechanisms to the radiative decay widths. All the units are in KeV. 




920 960 1000 1040 1080 920 960 1000 1040 1080 

Ms [MeV] Ms [MeV] 

(a) (b) 



Fig. 3. Left panel: fo — > p 7 decay width as a function of the scalar meson mass for the type-a loop. Solid line: n ir~ loop; 
dashed line: K + K~ loop. Right panel: ao — > uj"/ decay width as a function of the scalar meson mass for the type-b loop. Solid 
line: -K°r]Lup loop; dashed line: K + K~ K* + p° loop. 



dent on the scalar meson mass. This is specially relevant 
for the /o decay channels. 

In fig. [3] we show, as an example, the decay widths as a 
function of the scalar meson mass that one would obtain 
(without the convolution) for particular mechanisms for 

3(b)| mechanisms. 

de: 



the ty pe-a (fig. |3(a)[ ) and type-b (fig 
In fig 



one can see the strong dependence with the 
le type-a kaon loops close to the KK thresh- 



[3(a) 

energy of t 

old. This makes that small variations in the scalar meson 
mass change drastically the final radiative decay width. 
This is why it is very important to account properly for 
the scalar meson mass distribution by folding with the 
spectral function. For the type-a 7T7t loops and for the 
type-b, the dependence with the energy are smoother be- 
cause, as seen in table [3J they are dominated by 7T7t or 7777 
P1P2 mechanisms, which have the threshold far from the 
~ 980 MeV region of interest. Therefore, the influence of 
the convolution with the scalar mass distribution is less 
relevant for the type-b loops. 

At this point it is worth commenting about the rela- 
tive size of the divergent and convergent part of the type-b 
loops. We have checked that by using in the evaluation of 
the decay widths only the convergent (Eq. I25p or diver- 
gent (Eq.[26|) part, the contribution from the type-b loops 
would be from 3 to 10 times larger (depending on the 
channels) than those given in table [3J However, the sign 



of cZ^ v is opposite to that of cZ^ onv and then a strong 
destructive interference results leading to the type-b con- 
tributions shown in table [3] 

Finally, in table [4] we compare our result with other 
theoretical determinations. The model [I] (I) is based on 
quark models, ref. [I] (II) on kaon loops, ref. [23] (I) and (II) 
on VMD from direct SVP contact term and ref. [21] on 
meson loops by using the resonance Lagrangians of [3T] . 

We can see the wide dispersion among the different 
theoretical models. The previous calculations dealing with 
meson loops, refs. [1] (II) and [24], do not include the type- 
b loops considered in the present work. The model of 
ref. p](II) neglects the pion loops for the type-a mech- 
anisms which, as shown above and despite being small by 
themselves, affects the final result of the radiative decay 
widths in a significant manner after the coherent addi- 
tion to the dominant contributions. On the other hand, 
the model [T] (II) relies on an estimation of the gspp that 
should be viewed as qualitative [sic] [F and the model of 
[2~4"j uses couplings from the resonance Lagrangian of [3T] 
which we questioned before. On the contrary, the chiral 
unitary approach provides precise values for the couplings, 
including their phase, which have been successfully tested 
in many processes, and have uncertainties under control. 
In any case we should stress that all the results from the 
meson loops are of the same order of magnitude, and quite 
different from the quark and VMD models. 
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.1(1) 


.1(11) 


.23(1) 


23 (II) 


24 


present work 


fo - P U 1 


125 


3 


19 ±5 


3.3 ±2.0 


9.6 


4.2 ± 1.1 


fo -> W7 


14 


3 


126 ± 20 


88 ± 17 


15.0 


4.3 ± 1.3 


a -> p7 


14 


3 


3.0 ± 1.0 


3.0 ± 1.0 


9.1 


11 ±4 


ao — > luj 


125 


3 


641 ± 87 


641 ± 87 


8.7 


31 ± 13 



Table 4. Comparison of the radiative decay widths with other theoretical predictions. All the units are in KeV. 



The large dispersion among the theoretical models 
stresses the extreme sensitivity of these radiative decays 
to the theoretical models and on the nature of these res- 
onances. Hence, an experimental measurement would be 
very valuable to discern among theoretical models. The 
numerical values obtained here are within reach in present 
experimental facilities, in particular at COSY/ Juelich where 
an experiment is already approved [25] . 



4 Summary 

We have evaluated the radiative decays of the do (980) 
and /o(980) scalar mesons into a vector meson and one 
photon from the point of view of the chiral unitary ap- 
proach, where the ao(980) and /o(980) scalar mesons are 
dynamically generated by implementing unitarity in cou- 
pled channels in the pseudoscalar-pseudoscalar interac- 
tion. By evaluating the residues at the pole positions in 
unphysical Riemann sheets of the partial wave amplitudes, 
the couplings of the scalar mesons to the different PP 
channels can be obtained both in modulus and phase. 
Within this dynamical framework, the natural way to eval- 
uate the radiative decay widths is to consider the decay 
of the scalar mesons into their constituent PP building 
blocks and then allow the photon to be emitted from the 
pseudoscalar legs of the allowed loops. We consider not 
just the mechanisms with only pseudoscalar mesons in the 
loops, but also those with a vector meson in the loop. By 
using arguments of gauge invariance we show that the con- 
tribution of the loops involving only pseudoscalar mesons, 
which are in most of the cases the dominant ones, are 
convergent. The loops containing a vector meson are log- 
arithmically divergent, but this divergence can be recast 
into the two meson-loop function used in the scattering 
process, well under control, which leads to the dynamical 
generation of the scalar resonances. 

We make predictions for all the allowed /o/ao — * Vl 
decay widths including also an error analysis from the un- 
certainties in the parameters of the model. We show that 
even if some of the mechanisms are small by themselves, 
like the pion loops for the fo — ► p°7 or the loops contain- 
ing a vector meson in the loop, they affect strongly the 
final results due to non-trivial interferences with the dom- 
inant mechanisms. The sign of these interferences are well 
under control thanks to the knowledge of the phase of the 
couplings provided by the underlying unitary theory that 
generates dynamically the scalar resonances. 

These radiative decays are very sensitive to details 
of the models, hence an experimental measure would be 



highly valuable to discern among theoretical models. The 
COSY/WASA [26] scheduled experiment should be very 
useful in this respect. 
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d6c<iy 




A 




fo -» Pol 


7T + 7T~ 


-y/2/3 


-V2 




7T~7T~ + 


- v /2/3 


V2 






-1/V5 


-1/V2 




K K 


1 / /o 

-1/V2 


1/V 2 


fo -» 


K+K~ 




-W2 




K~K+ 


-1/^2 


1/V2 


Oo P07 


K+K~ 


-1/V2 


-1/V2 




K~K+ 


-1/^2 


1/V2 


ao -> W7 


K + K~ 


-1/V2 


-1/^2 




K~K+ 


-1/V2 


1/V2 




K+K° 


-1 


-1 



Table 5. Coefficients of Eq. (|14fl for type-a diagrams 



decay 


P1P2V2V1 


A 


5 


C 


/o -> Po7 




-y/2/3 


V2 




K + K-K*+p" 

UJ 


-1/V2 


1/V2 
1/V2 
1 


1/V2 


K-K+K*-p" 

UJ 


-1/V2 


1/V2 
1 


1/V2 


K»K°K* u p u 

UJ 


-1/V2 


-1/V2 

l/v^ 
1 


-1/V2 


K o K u K *» p u 

UJ 


-1/V2 


-1/V2 

1/V2 
1 


-1/V2 


/o -> ory 


-k + tt~ p + uj 


-y/2/3 


^2 




1T~ 7T + p~ UJ 


-y/2/S 


72 


V2 


n\°p°uj 


-a/2/3 


V2 


V2 


K + K-K"+p l) 

UJ 


-1/V2 


1/V2 
1/^2 
1 


1/V2 


K-K + K*-p" 

UJ 


-1/V2 


1/V2 
1/V2 
1 


1/V2 


K»K°K* u p u 

UJ 


-1/V2 


-1/V2 
1/V5 
1 


1/V2 


K°K u K* u p u 

UJ 

<5> 


-1/V2 


-1/V2 
1/V2 
1 


1/V2 



Table 6. Coefficients A, B and C of Eq. (|19[1 for the different allowed type-b diagrams. The continuation is in table 0. 
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decay 




A 


B 


C 


"o -> Pol 


K+K-K*+p u 

LU 


-1/V2 


1/V2 
1/V2 
1 


1/V2 


K~K + K*-p" 

LU 


-1/V2 


1/V2 
1/V2 
1 


1/V2 


K u K°K* u p u 

LU 


1/V2 


-1/V2 
1/V2 
1 


-1/V2 


K°K u K M 'p" 

LU 

<t> 


1/V2 


-1/V2 
1/V2 
1 


-1/V2 


■K° nP °LU 


1 


V2 


v/2/3 


r]TY°LULU 


1 


^2/3 


V2 


ao -> 


K+K-K'+p" 

LU 


-1/V2 


1/V2 
1/V2 
1 


1/V2 


K~K + K*-p" 

LU 


-1/V2 


1/V2 
1/V2 
1 


1/V2 


K u K°K* u p u 

LU 


1/V2 


-1/V2 
l/>/2 
1 


1/V2 


KOK^K^'p" 

LU 

4> 


1/V2 


-1/V2 
1/V2 
1 


1/V2 


Tv°riLup 


1 


^2 


^2/3 


ri7T°p°p 


1 


J ! 

v/2/3 


^2 




K+K°K*+p u 

LU 


-1 


1/V2 
1/V2 
1 


1 


K°K + K* u p u 

LU 


-1 


-1/V2 
1/V2 
1 


1 


ir + rip + u) 


1 




v/2/3 


r]1T + LULU 


1 


v/2/3 


72 



Table 7. ... continuation of table [6] 



